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Abstract 

The X(3872) seems to be a weakly-bound hadronic molecule whose constituents are two charm 
mesons. Its binding energy is much smaller than all the other energy scales in QCD. This separation 
of scales can be exploited through factorization formulas for production and decay rates of the X. 
In a low-energy effective field theory for the constituents of the X, the factorization formulas can 
be derived using the operator product expansion. The derivations are carried out explicitly for 
the simplest effective theory in which the constituents interact through a contact interaction that 
produces a large scattering length. The long-distance factors in the operator product expansions 
for various observables are calculated nonperturbatively in the interaction strength of the contact 
interaction. After renormalization of the coupling constant, all remaining ultraviolet divergences 
can be absorbed into the short-distance factors in the operator product expansions. 

PACS numbers: 12.38.-t, 12.39. St, 13.20.Gd, 14.40. Gx 
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I. INTRODUCTION 



The X(3872) is a narrow resonance near 3872 MeV discovered by the Belle collaboration 
in 2003 through its decay into J/ip'K^'K~ Jl|. Its existence was subsequently confirmed by 
the CDF, Babar, and DO collaborations [2,0j0|- Its mass Mx is extremely close to the 
threshold for the charm mesons D'^ and D*'^ [q]: 

Mx - (Mz)o + M^.o) = +0.7 ± 1.1 MeV. (1) 

The X(3872) is narrower than most of the known charmonium states 0]: 

Vx < 2.3 MeV (90% C.L.). (2) 

The observation of the decay X J/i'l implies that the X has charge conjugation quan- 
tum number C = + Pt]. Analyses of the discovery decay mode X J/-?/; vr+vr", including 
the angular correlations between the J/ip and the pions and the vr+vr^ invariant mass distri- 
bution, strongly favor the spin and parity quantum numbers = 1"*" These properties 
are compatible with the identification of X as a weakly-bound molecule whose constituents 

are a superposition of charm meson pairs SM Q H Q Q 111 H 113 El HH H H 

2l,|2l|23: 

X = ^ {D*'>D'> + . (3) 

v2 

If this identification is confirmed, the X(3872) would be the first unambiguously identified 
member of a new class of hadrons: mesonic molecules jH 27, H 29l. l30|. 



If the X(3872) is a weakly-bound mesonic molecule, it shares an important feature with 
the simplest baryonic molecule, the deuteron. Their binding energies are both small com- 
pared to the natural energy scale associated with the exchange of the lightest meson, the 
pion. That natural energy scale is m^/(2Mi2), where M12 is the reduced mass of the two 
constituents. The binding energy 2.2 MeV of the deuteron is small compared to the natural 
scale of about 20 MeV. The measurement of the mass of the X in Eq. (0) implies that its 
binding energy {M^o + M^i.o) - Mx is between -2.4 MeV and 1.2 MeV at the 90% con- 
fidence level. The small width in Eq. (j2I) further suggests that the mass of the X must be 
below the threshold for the charm mesons: Mx < M^io + M/j.o . Thus the binding energy of 
the X is small compared to the natural scale of about 10 MeV. The deuteron has an S-wave 
coupling to its constituents, the proton and the neutron. The quantum numbers J^^ = 
of the X implies that it also has an S-wave coupling to its constituents. The combination 
of the small binding energy compared to the natural energy scale and the S-wave coupling 
to the constituents implies that the deuteron and the X(3872) have universal properties 
that are determined by the large scattering length a of their constituents jl^J • The simplest 
example of a universal result is a simple formula for the binding energy of the molecule: 
Ex = l/(2Mi2a^). The universality of few-body systems with a large scattering length has 
many applications in atomic, nuclear, and particle physics 8]. The universal features of the 
X(3872) were first exploited by Voloshin to describe its decays into D^D'^ir^ and D'^D^'j, 
which can proceed through decay of the constituent D*^ or Z)*° [12]. Universality has also 
been applied to the production process T(4S') — > ti'^tt" + X [16], to the production process 



B ^ K + X |18|, |21|, to the line shape of the X [22], and to decays of X into J/i/j and pions 

3- 

The tiny binding energy of the X(3872) provides a new energy scale that is much smaller 
than the other scales in QCD, including the pion mass and the scale Aqcd associated 
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with nonperturbative effects. In Ref. [22], this separation of scales was exploited by using 
factorization formulas to separate certain observables into long-distance factors that involve 
only energy scales comparable to the binding energy and short-distance factors that involve 
all the higher energy scales of QCD. The long-distance factors can be calculated using an 
effective field theory for the constituents of the X that describes the lowest energy scale. 

In this paper, we show how the factorization formulas can be derived using the operator 
product expansion for the effective field theory that describes the constituents of the X. 
The effective field theory that describes the X(3872) is complicated by the spin 1 of the 
constituent D*^ and by its charge conjugation quantum number C = +, which implies 
that it is the superposition of D*^D^ and D^D*^ given in Eq. (jS)). Another complication 
is that the D^D^n^ threshold is only about 8 MeV below the D*^D^ threshold O. We 
will therefore illustrate the operator product expansion formalism using a simpler model in 
which these complications are absent. The simplest such model is a scalar meson model 
in which the constituents are spin-0 mesons with a contact interaction that gives a large 
positive scattering length a. The generalization to a realistic model with charm mesons is 
then straightforward. 

In Sec. m we define the minimal charm meson model that can describe the X(3872) and 
the simpler scalar meson model. In Sec. IIIH we explain how the operator product expansion 
can be used to separate scales in short- distance production and decay rates. In Sec. IIVI 
we give exact nonperturbative results for long-distance observables in the scalar meson 
model. In Sec. El we apply the operator product expansion to short-distance production 
and decay rates in the scalar meson model and we calculate the long-distance factors in the 
operator product expansion. We show that after renormalization of the coupling constant, 
all remaining dependence on the ultraviolet cutoff can be eliminated by renormalization 
of Wilson coefficients in the operator product expansion. In Sec. IVIl we show how the 
factorization formulas can be simplified by expanding in inverse powers of the large scattering 
length. In Sec. IVIIl we extend the results of Sees. IIVI IV| and EH to the charm meson model. 
We summarize our results in Sect. IVIIIl 

II. EFFECTIVE FIELD THEORIES 

In this section, we define the scalar meson model and the minimal charm meson model. 
Both of these effective field theories have an S-wave bound state. We introduce interpolating 
fields for the S-wave bound states. 

A. Scalar Meson Model 

We will illustrate the derivation of factorization formulas using the operator product 
expansion in a simpler model that we call the scalar meson model. This model has an S- 
wave bound state which we will refer to as X. The constituents of the bound state are scalar 
mesons Di and D2 with masses Mi and M2 satisfying Mi < M2. The scalar meson model is 
a nonrelativistic quantum field theory with two complex scalar fields Di{f,t) and D2{r,t). 
The free terms in the lagrangian are 

U„ = Dl (4 - M, + ^V=) D, + Dl (4 - + J^V^) . (4) 
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FIG. 1: Vertex for the D1D2 contact interaction. 

A superscript f on a field represents its complex conjugate. The interaction term in the 
lagrangian for the scalar meson model is 

Ant = -\oDIdId,D2 . (5) 

The coupling constant Aq has mass dimension —2. The subscript on Aq emphasizes that it 
is a bare coupling constant that depends on the ultraviolet cutoff on the momenta of the 
particles in loop diagrams. If the D\d\DiD2 interaction is treated nonperturbatively, there 
is an S-wave bound state. 

It is convenient to introduce concise notations for the reduced mass of the Di and D2 
and for the sum of their masses: 

Af., = Jli^. (6a) 

M1+M2 ^ ' 

Mi+2 = M1 + M2. (6b) 

If the scalar meson model is an effective field theory for a more fundamental theory in which 
the mesons Di and D2 interact by the exchange of other mesons, the natural momentum 
scale for low-energy processes is the mass m of the lightest meson that can be exchanged. 
If m <^ M12, the natural energy and momentum scales associated with the exchange of that 
meson are m? /M12 and m, respectively. We assume that the binding energy of the molecule 
X is small compared to the natural energy scale: 

Mi+2 - Mx < mVMi2 . (7) 

The scalar meson model describes the threshold region where the invariant mass M of Di 
and D2 is very close to M1+2: 

IM-M1+2I <mVMi2. (8) 

This constraint restricts the possible scattering states to D1D2. 

The bare coupling constant Aq of the scalar meson model must depend on the ultravi- 
olet cutoff A in such a way that low-energy observables are independent of A. There are 
many alternative renormalization prescriptions that can be used to eliminate the explicit 
dependence on A and Aq. One renormalization prescription is to eliminate Aq in favor of 
a renormalized coupling constant A. Another renormalization prescription is to eliminate 
Ao in favor of the scattering length a of the two heavy mesons. The scattering length can 
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be defined in terms of tlie T-matrix element for elastic D1D2 scattering with zero relative 
momentum: 

Ti2^i2 {p = 0) = -87rMi+2 a. (9) 

This is the T-matrix element for relativistically normalized particles in the initial and final 
states. If either Di or D2 is an unstable particle, there is an inelastic scattering channel for 
DiD2- This implies that the scattering length a has a negative imaginary part. The most 
convenient renormalization prescription for our purposes is to eliminate An in favor of the 
energy -Epoie at which the Green's function for D1D2 D1D2 has a pole 23|. That energy 
can be expressed in the form 

^poic = Mi+2-7V(2Mi2), (10) 
where 7 is the complex binding momentum: 

7 = 7re + ilim- (H) 

Unitarity requires 7im to be positive. We assume that 7re is also positive, in which case the 
energy Epo\e is the complex energy of the unstable bound state we denote by X. The real 
part of -Epoio defines the pole mass of the molecule: 

Mx = Mi+2 - ill - lL)/m,2). (12) 
The imaginary part of -Epoic multiplied by —2 can be interpreted as the width of the molecule: 

Tx = 27,e7im/Mi2. (13) 

The magnitude of the complex parameter 7 is assumed to be small compared to the natural 
momentum scale: I7I -C m. This implies the condition on the binding energy in Eq. ((7j). 



B. Minimal Charm Meson Model 

The charm mesons D° and D*^ with nonrelativistic energies and momenta can be de- 
scribed by a nonrelativistic quantum field theory with a complex spin-0 field D{f,t) and a 
3-component complex spin-1 field D{f,t). Their antiparticles and D*^ can be described 

by corresponding fields D{f,t) and I){r,t). The free terms in the lagrangian density for 
these particles are 

£„„ = (4 - M„. + ^V^) D + flt (4 - M„. + ^V^) D 

+5. . (4 _ M„.,. + _i_v^) S + flt . (4 - Mo. + ^V^) D .(14) 

The simplest interaction term that can produce an S-wave bound state in the C = + 
channel is 

£int = -Ao [dD + DD^ ^ ■ [dD + DD^ . (15) 

We will refer to the effective field theory with lagrangian given by Eqs. (fT^ and (fT^ as the 
minimal charm meson model. If the interaction in Eq. (P3|) is treated nonperturbatively, 
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there is a S-wave bound state with spin 1 that can be identified with the X(3872). The 
effects of decays of the X can be taken into account by allowing the coupling constant Aq in 
Eq. ()15|) to have an imaginary part. 

An ultraviolet cutoff A is required to regularize ultraviolet divergences generated by the 
interaction term in Eq. (jl5|) . The natural scale for the ultraviolet cutoff is the pion mass m^r. 
The bare coupling constant Aq must depend on A in such a way that low-energy observables 
are independent of the cutoff. There are many alternative renormalization prescriptions 
that can be used to eliminate the explicit dependence on A and Aq. For example, the 
complex parameter Aq can be eliminated in favor of a renormalized coupling constant A 
or in favor of the complex scattering length of the charm mesons. The most convenient 
renormalization prescription for our purposes is to eliminate Aq in favor of the mass and 
width of the X(3872) or equivalently the complex binding momentum 7. An alternative 
statement of this renormalization prescription is that the Green's function for D*^D'^D*^ 
has a pole at the energy -Epoie given by Eq. (fTI]|) . 



C. Interpolating fields for X 

In the scalar meson model, the local composite operator dIdI{x) has a nonzero amplitude 
to create X from the vacuum. Thus DiD2{x) can be used as an interpolating field for X. 
The resulting propagator for X is 

iAxiE,P) = y"rf^xe'^-"(0|A^2(x)/^I/^^(O)|0), (16) 

where P ■ x = P^x^ and P^ = {E,P) is the 4-momentum of the X. The propagator is 
a function of E and P = \P\. The Galilean invariance of the scalar meson model implies 
that it depends only on the combination E — P^/(2Mi2). Our renormalization prescription 
implies that this propagator at P = has a pole in E at the complex energy -Epoie given in 
Eq. (fTUI) . The behavior of the propagator near the pole defines a wavefunction normalization 
constant Zx'- 

^A^(i^,0)^— (17) 

Because the composite operator D1D2 has mass dimension 3, the propagator iAx{E, P) has 
mass dimension 2 and Zx has mass dimension 3. 

T-matrix elements involving X in the final state can be obtained from connected Green's 
functions involving the operator D1D2 by using the LSZ formalism \^^. The connected 
Green's function in momentum space with an external line associated with a -D1-D2 operator 
is amputated by multiplying by the inverse propagator for X, evaluated on the energy shell 
E = Epoie, and then multiplied by Z]/^ to obtain the T-matrix element. T-matrix elements 
involving X in the initial state can be obtained in a similar way from connected Green's 
functions involving the operator DIdI- 

In the charm meson model, X(3872) is identified as a bound state whose constituents are 
the C = + superposition of charm mesons in Eq. Q. A convenient interpolating field for 
the X is the local composite operator D^D{x) + DD^{x). The resulting propagator for the 
X is 

i^l{E, P) = j e'^■%^ [D'D{x) + DD\x)) {DW{0) + DD^(0))^ |0). (18) 
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FIG. 2: Vertices for the D1D2 and d\D2 operators. 



If P =0, the propagator has a pole aX E = -Epoie- Its behavior near the pole defines a 
normalization factor Zx'- 



zA^(E,0) 



iZ 



X 



E — Epoie + is 



■ 6'^ . 



(19) 



III. OPERATOR PRODUCT EXPANSION 

The scalar meson model defined by the lagrangian in Eqs. (@)) and (0) can be a low- 
energy approximation to a more fundamental Lorentz-invariant quantum field theory. If the 
fundamental quantum field theory includes a high-energy process that can create D1D2 with 
invariant mass near M1+2, that process will involve momenta ranging from the highest energy 
scale to momenta smaller than m. If long-distance effects involving momenta much smaller 
than m can be separated from short- distance effects involving momenta of order m and 
larger, we can use the scalar meson model to calculate the long-distance effects. Expressions 
for physical quantities in which short-distance effects and long-distance effects are separated 
into multiplicative factors are called factorization formulas. The tool required to separate 
long-distance effects from short-distance effects is the operator product expansion. 



A. Short-distance production processes 

The general production process for D1D2 has the form A B + D1D2, where A and B 
each represent one or more particles. There can be analogous production processes for X. 
We consider the production of -D1-D2 near their threshold, so that their invariant mass M 
satisfies the inequality in Eq. (jHJ. We also assume that the relative momentum p of the Di 
and D2 is small compared to m. We call A B + D1D2 a short- distance production process 
if all the particles in A and B have momenta in the -D1-D2 rest frame that are of order m or 
larger. This condition implies that the amplitude for A ^ B -\- D1D2 can be expanded in 
powers of the small energy difference M — Mi_|_2 divided by m^/Mi2 and larger energy scales 
and in powers of the small relative momentum p divided by m and larger momentum scales. 
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The operator product expansion can be used to express the T-matrix elements in the forms 



r[A^B + DiD2] = V4MiM2^Cf'"(DiA|a(a; = O)|0), 



(20a) 



n 



r[A^B + X] = v/2Mx^Cj'"(X|a(a; = O)|0). 



(20b) 



n 



The sums are over local operators On in the effective field theory. They can be restricted 
to operators with a nonzero matrix element between and the vacuum state |0). The 

arguments of the operator are the origin x = of space and time. The operator matrix 
elements are evaluated in the rest frame of D1D2 or X. We use the standard nonrelativistic 
normalizations for the states or (X| in the operator matrix elements. We use the 

standard relativistic normalizations for the initial and final states in the T-matrix elements. 
The factors of V4M1M2 and ^/2Mx account for the differences between the normalizations 
of the states in the operator matrix elements and the T-matrix elements. The Wilson 
coefficients C^'" in Eq. ()20|1 are functions of the 4-momenta and polarization 4-vectors of 
the particles in A and B and of the total 4- momentum P'^ of a D1D2 pair that is produced 
exactly at threshold with invariant mass Mi_|_2. They also depend on energy scales of order 
m^/Mi2 and higher and on momentum scales of order m and higher. The only dependence 
on whether the final state includes -D1-D2 or X is in the operator matrix elements. The 
leading terms in the expansions of the T-matrix elements in powers of M — M1+2 and p 
are the terms with the lowest dimension operator D|d|(0). In Feynman diagrams, the local 
operator D\d\ is represented by an open dot from which a Di line and a D2 line emerge, 
as illustrated in Fig. |21 The Feynman rule for this vertex is 1. 

The operator product expansions in Eqs. (^01) provide the desired separation of long- 
distance effects and short-distance effects only if the local operators On are chosen to be 
renormalized operators that have ultraviolet-finite matrix elements. The simplest local com- 
posite operators, such as D1D2, generally have matrix elements that are ultraviolet divergent. 
However the multiplicative renormalizability of local composite operators implies that the 
relation between the simple local operators On and their renormalized counterparts On 
can be expressed as 



where Z is an infinite-dimensional matrix of renormalization constants. The separation of 
short- distance effects and long-distance effects in the T-matrix elements in Eqs. (|20p is ac- 
complished by making the substitutions in Eq. (PT|) for the operators On- The long-distance 
factors are matrix elements of the renormalized local operators Om {x = 0). The short- 
distance factors are the coefficients of these matrix elements, which are sums of products of 
Wilson coefficients C^'" and renormalization constants {Z~^)nm.- We will find it more con- 
venient to work with simple composite operators rather than renormalized operators. The 
combination of the operator product expansion and the multiplicative renormalizability of 
these operators will be used to separate short- distance effects from long-distance effects. 

B. Short-distance decay processes 

The fundamental theory may also allow transitions -D1-D2 —>■ C from -D1-D2 with invariant 
mass Mc satisfying \Mc — M1+2I ^ w? /M\2 to a final state C that includes particles other 




(21) 



m 
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than Di and D2. If Mc = Mx, there can be analogous transitions X ^ C. If the sum of the 
masses of the particles in C is substantially smaller than M1+2, some of the particles in C 
must emerge with large momenta. We define a short- distance transition to be one for which 
all the particles in C have momenta in the center-of-momentum frame that are of order m 
or larger. The T-matrix element for such a process can be expanded in powers of the small 
energy difference Mq — M1+2 divided by m^/Mi2 and larger energy scales and in powers 
of the small relative momentum p of the Di and D2 divided by m and larger momentum 
scales. The operator product expansion can be used to express the T-matrix elements in 
the forms 

T[DiD2 ^C] = v/4MiM2^C^ (0|a(x = 0)\D,D2), (22a) 

n 

T[X^C] = V2Mx$^C^(0|a(x = O)|X). (22b) 

n 

The sums over local operators On of the -D1-D2 model can be restricted to those with a 
nonzero matrix element between the vacuum (0| and \D1D2). The only dependence on the 
initial states is in the operator matrix elements. The leading terms in the expansions of the 
T-matrix elements in powers of Mc — M1+2 and p are the terms with the lowest dimension 
operator DiD2{0). The complete separation of short-distance and long-distance effects is 
accomplished by using Eq. (j2ip to eliminate the operators On in favor of renormalized 
operators. 

C. Line shape 

If the fundamental theory includes short- distance processes that allow the production of 
X via A ^ B + X and the decay of X via X — ^ C, it also allows the process A ^ B + C, 
where C represents the same particles but with a variable invariant mass Mc instead of Mx- 
This process has a resonant enhancement when Mc is near the -D1-D2 threshold as specified 
by Eq. (jSJ. If each of the particles in A and B is well-separated in momentum space from 
each of the particles in C, the T-matrix element for this process can be described within the 
effective field theory by a double operator product expansion: 

T[A^B + C] = CY + ^Cf'"C^ I d^xe^^-^(0|O„(x)a(O)|0). (23) 

m,n 

The sum over operators Om can be restricted to those with a nonzero matrix element between 
(D1D2I and the vacuum state |0). The sum over operators 0„ can be restricted to those 
with a nonzero matrix element between the vacuum (0| and \DxD2}. The Wilson coefficients 
C^'"" and are the same ones that appear in the operator product expansions in Eqs. 
and ()22|). In the Fourier transform of the vacuum-to-vacuum matrix element in Eq. ()23p. 
the 4- vector is = {Mc,0). The leading term in the expansion of the T-matrix element 
in powers of Mc — M1+2 divided by m'^/M^ and larger energy scales is the term with the 
lowest dimension operators DiD2{x) and dIdIIO). The first term C^'*" on the right side of 
Eq. (j23|) takes into account the direct production of C at short distances. This term can 
be expanded in powers of the small energy difference Mc — M1+2 divided by m^/Mi2 and 
larger energy scales. The leading term in the expansion is a constant independent of Mc- 
The complete separation of short- distance and long-distance effects is accomplished by using 
Eq. ()21|) to eliminate the operators On in favor of renormalized operators. 
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XX-X-XIX 



FIG. 3: The amputated connected Green's function iAo{E) for D1D2 D1D2 at 0* order in g. 
It can be obtained by summing a geometric series of diagrams. 

IV. LONG-DISTANCE PROCESSES 

In this Section, we give the results for several quantities in the scalar meson model that 
depend only on long distances: the Green's function for -D1-D2 -D1-D2, the cross section 
for elastic -D1-D2 scattering, and the propagator for the bound state X. 



A. Green's function for D1D2 D1D2 

In the scalar meson model, all the observables for processes near the D1D2 threshold 
are related in a simple way to the Green's function for D1D2 — > D1D2. We denote the 
amputated connected Green's function for D1D2 — > -D1-D2 by iAq^E), because it depends 
only on the total energy E in the -D1-D2 rest frame. It can be calculated nonperturbatively 
by summing the geometric series represented by Fig. 01 to all orders in Aq: 

'■^»(^' - '''' 

where zLq^E) is the amplitude for the propagation of the -D1-D2 pair between successive 
contact interactions. The function Lq{E) has an ultraviolet divergence that can be isolated 
into an additive term that is independent of E. It can be expressed as 

Lo(^) = Lo(Mi+2) + ^V-2Mi2^i2-^£, (25) 
where E12 is the energy relative to the -D1-D2 threshold: 

E^2 = Mi+2. (26) 

The ultraviolet divergence is contained in the term Lo(Mi+2). If we use an ultravio- 
let momentum cutoff A ^ |Mi2-E'i2|^''^, this term has a linear ultraviolet divergence: 
Lo(Mi+2) = — (Mi2/7r^)A. If we use dimensional regularization, this term is Lo(Mi+2) = 0, 
because dimensional regularization sets power ultraviolet divergences equal to 0. 

One can identify the combination I/Aq — -Lo(Mi+2) on the right side of Eq. as a 
renormalized inverse coupling constant 1/A. One possible renormalization prescription is to 
eliminate I/Aq in favor of 1/A. Another possible renormalization prescription is to eliminate 
Ao in favor of the scattering length a defined in Eq. (jH)) . The most convenient renormalization 
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prescription for our purposes is to demand that the pole in the amphtude Ao{E) be at the 
energy -Epoie given in Eq. ()1U|). Inserting the expression for Lq{E) in Eq. ()25|) into the 
amphtude in Eq. and ehminating Aq in favor of 7, the amphtude reduces to 

ME) ^ (27) 

The renormahzed expression for the amphtude Aq{E) in Eq. ()27|) fohows from the renor- 
mahzation prescription in Eq. (|Tfl|) and is independent of the regularization scheme. The 
expression for the bare couphng constant is 

^ \ ('2g') 

° Lo(Mi+2)+Mi2(7rc + ^7im)/(27r)' ^ ' 

With dimensional regularization, Lo(Mi+2) = 0, so Eq. ()28|) gives a finite relation between 
Ao and the binding momentum 7: Aq = 27r/(Mi27). We will see later that that the naive 
use of dimensional regularization can be misleading. 



B. Elastic D1D2 Scattering 

We can use the amplitude in Eq. (j^7j) to determine the T-matrix element for the elastic 
scattering of Di and D2 with relative momentum p. In the D1D2 center-of-momentum 
frame, the total energy of the Di and D2 is 

E,M = Mi+2+pVi2Mu). (29) 

The energy variable E12 defined in Eq. (j26|) reduces to E12 = p^/(2Mi2). The T-matrix 
element is obtained by evaluating the amphtude Ao{E) in Eq. (|77j) at the energy Ecm{p) in 
Eq. ()29|) and multiplying by the factor 4M1M2 to account for the relativistic normalization 
of states: 

87rMi+2 , . 

Ti2^i2 P = -—■ (30) 

'J + tp 

The complex D1D2 scattering length defined by Eq. Q is therefore simply 

a = l/(7re + i7im)- (31) 

We obtain the cross section for elastic -D1-D2 scattering by squaring the T-matrix ele- 
ment, integrating over the phase space of the Di and D2 in the final state, and multiplying 
by a flux factor. The energy Ecm{p) in the D1D2 rest frame is assumed to be close to 
Mi+2, as specified by the condition in Eq. (jH)). The product of the phase space factor 
A^/2(^^^(p)^Mi,M2)/(87rM2^2) and the flux factor l/(4Mi+2p) can therefore be approxi- 
mated by l/(167rMf^2)- The cross section is 

The argument (p) of D1D2 in the initial state implies that the Di and D2 have momenta 
—p and +p, respectively. The momenta of the Di and D2 in the final state are not specified 
because they have been integrated over. 
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C. Propagator for X 



If the local composite operator DiD2{x) is used as an interpolating field for X, the 
propagator for X is given in Eq. ()16|) . The diagrams for the propagator of X are shown in 
Fig. El In the rest frame P = 0, these diagrams form a geometric series whose sum is 

This propagator has a pole in E at the same energy -Epoie given in Eq. ()10|) as the amplitude 
Ao{E) in Eq. ()24|) . Near the pole, the behavior of the propagator at P = is given in 
Eq. (fTTj) . Using Lo(ii^poie) = l/^o, we determine the wavefunction normalization factor to 
be 

^ Mi- <^^' 

Using the expression for Aq{E) in Eq. (j2D), the propagator for X in Eq. (jH^ can be 
expressed as 

Ax(i?,0) = -A(i?)^^. (35) 

The alternative expression for Aq{E) in Eq. (|27|) shows that, after renormalization of the 
coupling constant, it does not depend on the ultraviolet cutoff A. Thus the propagator in 
Eq. depends on A only through the factor Lq{E)/Xq. That there is some dependence 
on A is not a surprise, because we have used the simple composite operator DiD2{x) as the 
interpolating field for X rather than a renormalized operator. As we shall see in Sec. |3 
XoDiD2{x) is a renormalized operator whose matrix elements between the vacuum (0| and 
\D1D2) or \X) do not depend on the ultraviolet cutoff. The propagator for the renormalized 
operator XoDiD2{x) is obtained by multiplying the expression in Eq. (jH^ by Ag. But this 
propagator depends on A through the factor XqLq{E). To obtain a renormalized propagator 
that does not depend on A, one must add the A-dependent constant iXo to the propagator 
for the renormalized operator XqDiD2{x): 

iXl Ax{E, 0) + iXo = -iAo{E) . (36) 

Thus the renormalized propagator is essentially just the amplitude Ao{E) in Eq. (j27p . The 
need for adding the constant term in Eq. (j36|) is related to the fact that if an external source 
coupled to a composite operator is added to the lagrangian, renormalization sometimes re- 
quires the addition of terms with higher powers of the source js^]- For example, the addition 
of the term DiD2 + h..c. creates new ultraviolet divergences that can only be cancelled by a 
J^J term. Such a term is required even in the absence of any interactions. To implement the 
LSZ prescription for T-matrix elements for processes with X in the initial or final state, it is 
not necessary to use a renormalized propagator. We will use the unrenormalized propagator 
for X in Eq. (jSH)) for this purpose. 



V. SHORT-DISTANCE PROCESSES 

In this Section, we consider processes in the scalar meson model that involve both short- 
distance and long-distance effects: short-distance production rates, short-distance decay 
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FIG. 4: Diagrams at 0* order in g for the vacuum-to-DiD2 matrix element of the operator DjDg- 

rates, and the line shape of the bound state in a short- distance decay mode. We use the 
operator product expansion to derive factorization formulas in which those short- distance 
effects and long-distance effects are separated. After renormalization of the coupling con- 
stant, all remaining depedence on the ultraviolet cutoff can be removed by renormalization 
of the Wilson coefficients in the operator product expansion. 

A. Short-distance Production of X and D1D2 

We consider the short-distance production processes A ^ B + D1D2 and A ^ B + X, 
where A and B both represent one or more particles whose momenta in the D1D2 or X 
rest frame are all of order m or larger. The operator product expansions of the T-matrix 
element for such short-distance processes are given in Eqs. ()20|) . The leading terms in the 
expansions are the ones with the local operator D\d\{^). If we keep only these leading 
terms, the Lorentz-invariant T-matrix elements reduce to 



We first consider the T-matrix element for the production of D2D2 via the short- distance 
process A B+D1D2. We take the Di and D2 in Eq. ()37a|l to have relative momentum p in 
the -D1-D2 rest frame and unspecified total momentum. Their invariant mass is M = Ecm{p), 
where i?cm(p) is given in Eq. This invariant mass is assumed to be close to M1+2, 

as specified by the inequality in Eq. (jH)). The Feynman diagrams for the vacuum-to--Dii^2 
matrix element, which are shown in Fig. |3J form a geometric series. The matrix element in 
Eq. ()37a|) is therefore 



The argument (p) of the -D1-D2 state implies that the Di and D2 have momenta —p and 
+p, respectively. The sum of the diagrams in Fig. H] differs from the geometric series of 
diagrams for iAo{E) in Fig. El only by the multiphcative factor l/(— iAo). Since Ao{E), 
which after renormalization of the coupling constant is given by Eq. ()27|) . does not depend 
on the ultraviolet cutoff, the operator XoDIdIIx) is a renormalized operator whose matrix 
elements do not depend on the ultraviolet cutoff. The T-matrix element in Eq. (j37ap can be 
separated into a short-distance factor and a long-distance factor by taking the long-distance 
factor to be the matrix element of the renormalized operator: 



r[A^B + D^D2] = y^AM^2C^'^^ {DiD2\dIdI{0)\^) 



(37a) 
(37b) 



T[A^B + X] = v/2MxCf'''(X|DjD|(O)|0). 




(38) 



T[A^B + DMp)] = - v^iMiMs (CJ^VAq) ^o(^cm(p)). 



(39) 
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FIG. 5: Feynman diagrams for the X propagator at 0* order in g. The interpolating field for the 
X is DiD2{x). 



After renormalization of the coupling constant, Ao{Ecra{p)) is given by the expression in 
Eq. (j2Z|), which does not depend on the ultraviolet cutoff A. The T-matrix element in 
Eq. will not depend on A if the short- distance factor C^'^^/ Aq does not depend on A. 
Equivalently, the dependence on A can be removed by a multiplicative renormalization of 
the Wilson coefficient C^'^^. 

We next consider the T-matrix element for the production of X via the short- distance 
process A B + X . The operator product expansion of the T-matrix element for this 
process is given in Eq. ()37b|) . The vacuum-to-X matrix element can be obtained via the 
LSZ formalism from the connected Green's function for the operator D\D2 that acts on the 
vacuum and an operator -D1-D2 associated with the X in the final state. This connected 
Green's function is identical to the propagator iAx(-E',0) given in Eq. (jHSj)- The matrix 
element in Eq. (j37b|l is the normalized on-shell amputated connected Green's function. The 
Green's function is amputated by multiplying by the inverse propagator [i/S.x{E, 0)]~^, which 
simply gives 1. In the rest frame of X where its momentum is P =0, the Green's function 
is put on shell by setting the energy E equal to the energy -Epoie in Eq. (fTUI) . although the 
absence of any dependence on E makes this condition moot. Finally the Green's function 

1 /2 ____ 

is normalized by multiplying by the factor , where Zx is given in Eq. ()34|) . Thus the 
vacuum-to-X matrix element is simply 

(X|DlDt(o)|0) = i^p:)^. (40) 

The only factors in the T-matrix element in Eq. ()37b|l that are sensitive to short distances 
are the Wilson coefficient C^'^^ and the factor of I/Aq from the matrix element. The T- 
matrix element in Eq. ()37b|) can be expressed as the product of a short- distance factor and 
a long-distance factor: 

T[A^B + X] = {C^'/Xo) (41) 

The short-distance factor C^'^^/Aq is the same one that appears in Eq. (jH^ . 

The factored expressions for the T-matrix elements in Eqs. (jH^ and ()41|) lead to factored 
expressions for the production rates. The rates for producing X and -D1-D2 are obtained by 
squaring the amplitudes and integrating over the appropriate phase space. If A consists of 
a single particle, the decay rate into B + X can be expressed in the factored form 

nA^B + X]=T^^^\^l + ^l\'/\ (42) 
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where is a short-distance factor with dimensions of mass: 



(43) 



The 4-momentum Pi_^_2 is that of a D1D2 pair exactly at threshold with invariant mass M1+2. 
We have chosen the long-distance factor in Eq. ()42|) to be the square of the long-distance 
factor in the T-matrix element in Eq. ()4H1 multiplied by M12 to make it dimensionless. All 
factors from integrating over the phase space of the particles in the final state are included 
in the short-distance factor. The differential rate for producing D1D2 with invariant mass 
M = Ecniip) given by Eq. (jSHl) can be expressed in the factored form 

[A^B + D.D^m = , .^^ , (44) 



dM lire + iilim + P)\ 

The short-distance factor is the same as in Eq. (021) • The long-distance factor is the 
product of 4MiM2|v4.o(-Ecm(p))P from the T-matrix element in Eq. the phase space 

factor \^/'^{Ecraip), Ml, M2) /(SnEcmipY), the kinematic factor Ecm{p)/'^ associated with 
the differential dM, and a factor Mi2/(2Mi+2) to compensate for the choice of prefactor 
in Eq. ()43p. Using the expression for the invariant mass Ecm{p) in Eq- ()29|) and replacing 
Ecraip) by Mi+2 in every factor that is insensitive to p, the product of the phase space and 
kinematic factors can be reduced to p/(47r^). 



B. Short-distance Decay of X 

The fundamental theory may have short-distance processes though which X can decay. 
We consider the short-distance decay process X — > C, where C represents two or more 
particles whose momenta in the X rest frame are of order m or larger. The operator product 
expansion for this decay process is given in Eq. ()22b|) . The leading term in the expansion is 
the one with the operator DiD2{x). If we keep only this term, the T-matrix element reduces 
to 

r[X -^C] = y^M^Cg (01^1^2(0) (45) 

The X-to-vacuum matrix element can be calculated in a similar way to the vacuum-to-X 
matrix element in Eq. ()4Up :^ 

f27r^)i/2 
M12A0 

This matrix element depends on the ultraviolet cutoff only through the factor I/Xq. In 
the T-matrix element in Eq. the separation of short-distance and long-distance effects 



^ Our notation might suggest that the matrix element in Eq. (|46|l is the complex conjugate of the matrix 
element in Eq. (|40|l . However \X) in Eq. (|46|l is an in state, while {X\ in Eq. (|40|) is the hermitian 
conjugate of an out state. These two states are related by the S-matrix: \X, out) ~ S\X, in). 
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can be accomplished by taking the long-distance factor to be the matrix element of the 
renormalized operator Ao-Di-D2(0): 



r\X ^C\ = v/2Mx (Cf,/Ao) (47) 



This T-matrix element will not depend on the ultraviolet cutoff A if the short- distance 
factor Cj^/Ao does not depend on A. Equivalently, the dependence on A can be removed by 
a multiplicative renormalization of the Wilson coefficient C^. 

The decay rate of X into the particles represented by C is obtained by squaring the T- 
matrix element and integrating over the phase space of those particles. It can be expressed 
in the factored form 

r[x-c] = r^^KU7,ir/\ (48) 

where F"-^ is a short- distance factor with dimensions of mass: 

- h / n (4% - 

We have chosen the long-distance factor to be the square of the long-distance factor in the 
T-matrix element in Eq. (j47p multiplied by M12 to make it dimensionless. 

The separation of short-distance and long-distance effects for the transition 
Di{~p)D2{p) — >■ C can be accomplished in a similar way. The operator product expansion 
for the T-matrix element is given in Eq. ()22a|) . The T-matrix element can be expressed as 
the product of the same short-distance factor as in Eq. (P7|) and a long-distance factor that 
includes a factor of Ao{Ecm{p)) , where Ecm{p) is the energy in Eq. ((23). 



C. Line shape of X in a Short-distance Decay Mode 

If the fundamental theory includes short- distance processes that allow the production of 
X via A ^ B + X and the decay of X via X ^ C, it also allows the process A ^ B + C , 
where C represents the same particles but with a variable invariant mass Mc instead of Mx- 
We assume that Mc is near the D1D2 threshold as specified by Eq. (jH)). If every particle in 
A and B has momentum in the rest frame of C of order m or larger, the T-matrix element 
for this process can be expressed as the double operator product expansion in Eq. (j^ . The 
Wilson coefficient C^'*" can be expanded in powers of Mc — M1+2 divided by m? /M12 and 
higher energy scales. The leading term in this expansion is simply a constant independent 
of Mc- The leading terms in the double sum come from the operators O„(0) = D\d\{{]) 
and Om{x) = DiD2{x). According to Eq. (fT^ . the Fourier transform of the matrix element 
of these operators in the C rest frame is just the X propagator, which is given in Eq. ()33|) 
or (jnSl), evaluated at E = Mc- Thus the T-matrix element reduces to 

n^^BH-Cl^C-.crcg^^^^. (50) 

The Wilson coefficients and the factor Lq{E)/Xq in the X propagator depend on the ultra- 
violet cutoff A. All the dependence on the energy can be isolated in a term that does not 
depend on A by using the fact that the combination in Eq. does not depend on A. 
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By subtracting and adding zC^'^^CJ^/Aq to the two terms on the right side of Eq. (j^Ujl . the 
T-matrix element can be expressed as 



c). (51) 

After renormahzation of the couphng constant, Ao{E) is given in Eq. (j27|) . The short- 
distance factors C^'^^/Ao and C^/Aq in Eq. cannot depend on A, because otherwise 
the T-matrix elements in Eqs. (jHIJj) . (PT|) . and (jTTj) would depend on A. Thus the T-matrix 
element in Eq. will not depend on A if the constant term C^'^ — iC^'^^Cj^/Ao does not 
depend on A. Equivalently, the dependence on A can be removed by an additive renormal- 
ization of the Wilson coefficient C^'*". It is convenient to express the T-matrix element in 
Eq. (|5ip in the form 



B,12 



/Ao) (Cfa/Ao) MMc) - {2'n/M^2)c 



B,C 
A 



(52) 



where complex constant with dimension of length that is completely determined 

by short-distance factors. The natural scale for c^' is 1/m, where m is the mass of the 
lightest meson that can be exchanged between Di and D2. 

The factored expression for the T-matrix element in Eq. (jH^ leads to a factored expression 
for the rate for A ^ B + C . The invariant mass distribution of the particles in C is obtained 
by squaring the T-matrix element and integrating over the momenta of all the particles 
in the final state. It is convenient to express the phase space integral in an iterated form 
corresponding to the production of the particles in B and an effective particle of mass Mc 
followed by the decay of that effective particle into the particles in C. If A is a single particle, 
the decay rate is 



T[A^ B + C] 



1 



2Ma 



n 



27r 
(27r)32E, 



C 



(27r)32E, 



c 



n 

ieB 



-^(2vr)^5W(Pj^-P^-Epr) 



■(27r)V^)(P, 



\T[A^B + C][ 



(53) 



The invariant mass Mq of the particles in C can be replaced by M1+2 everywhere except 
in the long-distance factor of the T-matrix element. In that long-distance factor, it can be 
expressed as 



Mc = Mi+2+pI/{2Mu) 



(54) 



where can be positive or negative. The variable pc is pure imaginary if Mq < M1+2 and 
real and positive if Mc > Mi_|.2- The differential decay rate of the particle A for Mc near 
Mi+2 reduces to 



dr 



dM, 



c 



[A^B + c] = (rfr^)2 



vr 



The short-distance factors and F*^ 



B,C 



ilre- \Pc\)+hi> 
1 

7re + i{lim+Pc) 



+ C 



are the same as in Eqs. 



other short-distance effects are contained in the complex constant 
distribution in Eq. is continuous at Mc = Mi+2- 



B,C 



Mc < Mi+2, (55a) 
Mc > Mi+2. (55b) 



and 



All 



The invariant mass 
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The separation of scales represented by the renormahzed operator product expansions 
for the T- matrix elements in Eqs. (HH), (jlZj), and can be obscured by using 

dimensional regularization. In a generic regularization scheme with ultraviolet cutoff A, 
the short-distance quantities Aq, C^'^^, Cf^, and C^'*" are insensitive to 7. They depend 
on A in such a way that the combinations C^'^^/Xq, C^/Aq, and C^'*" — iC^'^^C^/Ao do 
not depend on A. Dimensional regularization sets power ultraviolet divergences to zero. In 
particular, it sets Lo(Mi+2) = 0, so the expression for the bare coupling constant in Eq. (pHj) 
reduces to Aq = 27r/(Mi27). The Wilson coefficients C^'^^, Cj^, and C^'*" do not depend on 
the ultraviolet cutoff of dimensional regularization. Compatibility with other regularization 
schemes requires however that they depend on 7 in such a way that the combinations 
C^'^'^/Xo, C12/X0, and C^'*" — iC^'^^CjJ^/Ao are insensitive to 7. This requires that C^'^^ 
and Cf'2 have multiplicative factors of 7" and that C^' have an additive term with a 
factor 7~^. Thus the Wilson coefficients in dimensional regularization are not short-distance 
factors. Their dependence on 7 is similar to the dependence of the Wilson coefficients on 
A in other regularization schemes. That dependence cancels in the combinations of Wilson 
coefficients and Aq that appear in the renormalized operator product expansions for the 
T-matrix elements. 



VI. LARGE SCATTERING LENGTH EXPANSION 

Effective field theories can exploit a large separation of momentum scales by providing a 
simpler description of the lowest momentum scale. Another important feature of effective 
field theories is that they provide a systematic framework for improving the accuracy of 
the description to any desired order in the ratio of the small momentum scale and higher 
momentum scales. In this Section, we discuss how the accuracy of the results for the scalar 
meson model in Sees. IIVI and |3 can be systematically improved using an expansion in the 
large scattering length. We also explain how that expansion can be exploited to simplify 
some of the results in Sec. El 

In the scalar meson model, the smallest momentum scale is the scale I7I associated 
with the large scattering length. In a more fundamental theory, there may be many larger 
momentum scales, but the most important at low energies is the mass m of the lightest 
meson that can be exchanged between Di and D2. The model defined by the lagrangian 
in Eqs. (jl)) and (0) reproduces all effects that are not suppressed by powers of |7|/m. The 
model can be systematically improved so that it reproduces all corrections to any desired 
order in \'y\/m. We will discuss only the improvements required to reproduce corrections 
through first order in \^\/m. 

The only improvement in the effective theory that is required to decrease the errors to 
second order in |7|/m is to take into account the effective range for S-wave scattering. This 
parameter can be defined by the low-momentum expansion for the inverse of the T-matrix 
element: 

1 1 A . 1 2 \ 

T{p) 87rAfi+2 V*^ 2 / 

If we impose the renormalization condition that the Green's function for D1D2 — > D1D2 
has a pole at the energy -Epoie given in Eq. (fTUj). one expression that will give the correct 
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effective range is 



where p = iyj —2M12E12 — is. The corresponding T-matrix element for elastic -D1-D2 scat- 
tering is then 

In short-distance observables, there are additional terms in the expansions in |7|/m com- 
ing from higher dimension operators in the operator product expansion. For each addi- 
tional gradient in the operator, the operator matrix element will have an additional fac- 
tor of order 7. The dimensions from these additional factors of 7 must be compensated 
by factors of 1/m in the short- distance coefficients. Only operators with a single gradi- 
ent can give contributions that are suppressed by one power of |7|/wi. The matrix ele- 
ment of the operator V*(-Di-D2) between \DiD2{'p)) or \X) and the vacuum (0| vanishes in 
the center-of-momentum frame. The other independent operator with a single gradient is 
V'DiD2 - DiV'D2. The matrix element {^V'DlD2 - DiV'D2\X) must vanish because the 
operator is a vector and there are no vectors associated with the state \X) in its center-of- 
mass frame. The matrix element (0|V*-Di-D2 ~ Di'V^D2\DiD2{p)) is nonzero and propor- 
tional to p*. This operator gives a term in the T-matrix element for A — >■ B + DiD2{p) 
in Eq. that is linear in the momentum p but has a suppression factor of 1/m in the 
short- distance coefficient. Higher dimension operators in the operator product expansion 
will contribute to the T-matrix elements for A ^ B + X in Eq. ()41|). for X C in Eq. ()47p. 
and for A ^ B + C in Eq. ()52j) only at second and higher orders in \'y\/m. 

The systematic expansion in powers of |7|/m can be used to simplify the leading order 
results for short- distance observables. The T-matrix element for A ^ B + C in Eq. (I52p 
has a resonant term Ao{Mc) and a nonresonant term c^' . The resonant term Ao{Mc) 
includes a factor that is of order I/I7I when pc is of order |7|. The nonresonant term c^'*" 
is completely determined by short-distance effects, so the natural scale for c^' is 1/m. For 
Pc of order I7I, this amplitude is suppressed by |7|/m compared to the resonant term in 
Eq. (j52j) . One can therefore set c^'*" = by truncating the expansion at leading order in 
|7|/m. The invariant mass distribution in Eq. (|55|) then reduces to 

[A^B + C] = (r^F^) , ^"".^ , , Mc < Mi+2, (59a) 
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7re+ (7im+Pc)^ 



im 



i^'n .,2^u\^^^2 Mc > M,+2. (59b) 



This simple factorization formula was first derived in Ref. |22]. If the nonresonant amplitude 
c^' in Eq. (joojl is included, the systematic expansion in powers of |7|/m requires that all 
other terms that are first order in |7|/m also be included. This requires that the effective 
field theory be improved so that it takes into account the effective range. 

The above derivation of the simple factorization formula in Eq. ()59|) is much cleaner than 
the derivation in Ref. 2^1 . In Ref. , the authors used an ultraviolet momentum cutoff A. 



They obtained results that did not depend on A by taking the limit A ^ 00. The expression 
for the bare coupling constant Aq in Eq. (PH|) has a term Lo(Mi+2) in the denominator. Since 
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Lq{Mc) and Lo(Mi+2) are both linear in A, the product XoLq{Mc) approaches 1 in the hmit 
A — ^> oo. In this hmit, the last factor in the second term on the right side of Eq. (jSUj) reduces 
to 

^^^^^^-^MMc). (60) 



1 - XoLoiMc) Ao 

The factors of I/Aq can be combined with the Wilson coefficients C^'^^ and to obtain 
short- distance factors with finite limits as A ^ oo. In Ref. jl^l, the authors omitted the 
C^' term in Eq. (jSUj) . This gave the simple factorization formula in Eq. ()59p . In retrospect, 
omitting the C^' term in Eq. (pOj) can be justified by the observation that the natural scale 
for the coefficient c^'*" in Eq. is l/m-,,. If is identified with the ultraviolet cutoff 
A, then c^'*" — > in the limit A — > oo. The derivation in Ref. jHI blurred the distinction 
between the arbitary unphysical ultraviolet cutoff A, which can be taken to oo, and the 
physical short- distance scale m^, which is fixed. By maintaining the distinction between A 
and m^, we were able to separate the renormalization of the operator product expansion 
from the expansion in inverse powers of the large scattering length and give a much cleaner 
derivation of the factorization formula. 



VII. MINIMAL CHARM MESON MODEL 

In this section, we generalize the results of Sees. IIV| IV| and I VII for the scalar meson model 
to the minimal charm meson model defined by the lagrangian in Eqs. (fT^ and The 
bound state in this model is identified as the X(3872). For simplicity of notation, we denote 
the masses of -D° and D*^ by Mi = M^o and M2 = M^i.o. Thus M12 is the reduced mass of 
D° and D*'^ and Mi_|_2 is the sum of their masses. 



A. Long-distance processes 

The amplitude Lq{E) in the charm meson model for the propagation of D*^D^ or D^D*^ 
between contact interactions is given by the same expression in Eq. ()25p as in the scalar 
meson model. The Green's function for D*^D^ D*^D^ is diagonal in the vector indices of 
the spin-1 mesons. The diagonal entries are 

This differs from the expression for iAo{E) in Eq. ()24|) only in the factor of 2 multiplying 
Lq{E), which accounts for the fact that the particles in each of the loops in Fig. El can be 
either D*^D^ or D^D*^. The Green's functions for D*^D^ D^D*'^, D^D*^ D*^D°, 
and D^D*^ D^D*^ are also given by this same expression. If we use the renormalization 
prescription that the Green's function in Eq. (jHT|) has a pole in E at the energy -Epoie given 
in Eq. (fTn|l . the expression for the diagonal entries of the Green's function can be reduced 
to 

AiE) = (62) 

-7 + V-2Mi2^i2 - ie ^ ' 
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The complex parameter 7 determines the mass and width of a bound state with spin 1 that 
we identify as the X(3872). 

The Green's function A{E) in Eq. (jU^ differs from Ao{E) in Eq. ()27|) by a factor of 1/2. 
The T-matrix element for D*^D^ D*^D^ therefore differs from the expression in Eq. (|3Up 
by a factor of 1/2. The resulting expression for the cross section for elastic D*^D^ scattering 
therefore differs by a factor of 1/4 from the cross section in Eq. (j32j) for the charm meson 
model: 

^ (63) 

The argument (p ) of D*^D^ in the initial state implies that the D*^ and have momenta 
+p and —p, respectively. In the final state, the relative momentum of the D*^ and have 
been integrated over. The cross section for elastic D^D*^ scattering and the cross sections 
for D*^D^ D^D*'^ and D^D*'^ —>■ D*'^D'^ are also given by the expression on the right side 
of Eq. dni. 

If the local composite operator D^D{x) + DD\x) is used as the interpolating field for the 
X, the propagator for X is given in Eq. ()18|). The Feynman diagrams for the propagator 
iA^^{E, 0) in the charm meson model differ from those in Fig. El for the scalar meson model 
only in that each loop receives contributions from two pairs of particles, D*^D^ and D^D*^. 
Thus the diagonal entries of the X propagator can be obtained from the propagator in 
Eq. dnni) by replacing Lo{E) by 2Lo{E): 



The normalization factor defined by Eq. (fT^ is 

7r7 



(65) 



This differs by a factor of 2 from the normalization factor Zx in Eq. (j34p for the scalar 
meson model. 

The minimal charm meson model is an effective field theory that exploits the small ratio 
between the scale I7I associated with a large scattering length and all the larger momentum 
scales of QCD. At very low energy, the most important of these larger momentum scales is 
the mass m^r of the pion. The minimal charm meson model takes into account all effects 
that are not suppressed by powers of \'y\/m.,^. The model can be systematically improved 
so that it incorporates corrections to any desired order in \'y\/mT^. We will discuss only the 
improvements required to take into account corrections that are first order in |7|/m^. 

At first order \'~f\/mTj, it is necessary to take into account not only the large scattering 
length a+ in the C = + channel but also the effective range r+. It is also necessary to 
take into account the scattering length a_ in the C = — channel. These parameters can 
be defined by low-momentum expansions of the T-matrix elements analogous to Eq. 
Since r+ and a_ both have dimensions of length, the natural scales for these parameters are 
order 1/m^. If we impose the renormalization condition that the C = + channel amplitude 
has a pole in the energy at -Epoie given in Eq. (jlOj) . the Green's functions in the two channels 
can be written as 

^-'^> = h + .rt[l-r,h->rt/21 • (^«^' 

A-{E) = -(27r/Mi2)a_ , (66b) 
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where p = i\/ —2M12E12 — is in Eq. ()66a|) . The expression for the cross section for elastic 
D*0£)0 scattering that replaces Eq. (jUHjl is 



(7 + ip)[l -r+(7-ip)/2] 

The cross section for elastic D^D*^ scattering is given by the same expression. The cross 
sections for D*^D^ D^D*^ and D^D*^ D*^D^ are given by the same expression except 
that a_ is replaced by — a_. 



B. Short-distance production of X, D*^D^, and D^D*^ 

The operator product expansion for the charm meson model can be used to separate 
the rates for short-distance production and decay processes into short-distance factors and 
long-distance factors. We first consider the short-distance production process A B + X 
and the corresponding production processes for D*^D^ and D^D*^. A specific example of 
such a process is the discovery production process — > +X. The leading terms in the 
operator product expansions for these processes are those with the operators D*^Z)^(0) and 
D^l)*^(0). The expressions for the T- matrix elements analogous to Eqs. fj37|) are 



+ Cf' (D*°D°|Z}t^'t(o)|0)) , (68a) 
r[A^ B + D^D*^] = V4M1M2 (C^'' {D^D*°\D'^D^{0)\(D) 

+ C^'' (D°D*°|D^^'^(O)|0)) , (68b) 
T[A^B + X] = ^/2I^ (C^'' {X\D'W\0)\(lS) + C^'' {X\DW'\0)\(lS)y (68c) 
The matrix elements between the vacuum and the charm meson states are 

(D-»D°(p,m)|Z)ti)-'(O)|0) = ^;^°f ;;<P'> e-'(m), (69b) 

{D'D*'{p,mWW\om = ^f'^f-}:^^] e"(m), (69c) 

1 - 2\qLq[Ec^[p)) 

{D'D*%p,m)\DW^\0m = (l + Y^|^^§^^) ^"M, (69d) 

where e*(m) is the polarization vector for the spin-1 meson. The arguments {p,m) of the 
D*^D° and D^D*° states imply that the spin 1 and spin mesons have momenta +p and 
—p, respectively, and that the spin 1 meson has spin quantum number m. The matrix 
elements between the vacuum and the X are 

(X(m)|D^t^t(o)|0) = lz]l^e*'{m), (70a) 

{X{m)\DW\0)\(lS) = ^Zl/^e*\m), (70b) 
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where e*(m) is the polarization vector for the X and the normahzation constant Zx is given 
in Eq. ()65p . The factors of 1/2 in Eqs. (fTOj) come from the fact that the Green's function 
for the operators D^D{x) + DD'^{x) and D*^Z)"''(0) or D^i)*"'"(0) are equal to the propagator 
for the X given in Eq. (j64j) multiplied by 1/2. The T-matrix elements in Eqs. (|69p can be 
expressed in a form in which short- distance and long-distance effects are separated: 



'AM1M2 



2An 



4 I ^ A 



2An 



T[A^5 + X(m)] = v^2M 



X 



^A ~ ^A 



BA 



vl/2 



2Ao 



Ml 



12 



(m). 



(71a) 



(71b) 
(71c) 



-,B,: 



+ Cf'*)/Ao and 



These T-matrix elements do not depend on the ultraviolet cutoff if (C^ -r 
C^'* — C^'* do not depend on A. Equivalently, their dependence on A can be eliminated by 
renormalizations of the Wilson coefficients C^'* and C^'*. 

Since the T-matrix element for A ^ B + X in Eq. ()71cp is the product of a short- 
distance factor and a long-distance factor proportional to 7^''^, the rate can be expressed as 
the product of a short-distance factor and a long-distance factor proportional to |7|. If A 
consists of a single particle, the decay rate is 



T[A-^ B + X] = T 



B 



27V 



M 



I 2 I 2 

I7re + 7 



1 1/2 



12 



(72) 



We have chosen the long-distance factor to be the same as in Eq. 

The expressions for the invariant mass distributions for D*°D^ and L>°L>*° that follow 
from the T-matrix elements in Eqs. (|71ap and ()71b|) are much more complicated and depend 
on the types of particles in B. However these expressions simplify if we keep only the leading 
terms in the expansions in \'y\/m. For p of order I7I, the resonant amplitude A{Ecra{p)) in 
Eqs. fl71a|) and ()71b|) has a factor of order 1/|7|. The nonresonant terms in Eqs. ()71ap 
and ()71b|) involve only short-distance factors and are insensitive to the scale |7|. They are 
therefore suppressed relative to the resonant terms by \'~f\/m. If we take p to be of order I7I 
and keep only the leading terms in |7|/m, the invariant mass distributions reduce to 



dT 
dM 
dT 
dM 



[A^B + D*^D%p)] = r 



^B 

A 72 



p 



7re + (7im+p)^ 



[A-^ B + D''D*'\p)] = V 



^B 

A 0,2 

re 



P 



7re + (7im+P)^ 



(73a) 
(73b) 



We have replaced Ecm{p) by M1+2 everywhere except in the long-distance factor. The 
short- distance factors are the same as in Eq. ()72|1 . They cancel out of the ratio between 
Eq. ()73a|) or Eq. ()73b|) and Eq. ()72j) . This ratio differs from the ratio between Eq. and 
Eq. ()42|) in the scalar meson model by the probability 
the C = -|- channel. 



i for the and D^D*^ to be in 
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The factorization formulas in Eqs. ()72p and ()73p were first derived in Refs. [l8[ and 21 



for the case 7im = 0. They were apphed to the decays of B mesons to K + X, K + D*^D^, 
and K + D^D*^ . The factorization formulas were generalized to the case 7im > in Ref. 2^ |. 



The short- distance coefficient in Eqs. (j72|) can be eliminated in favor of the decay rate 
T[A^ B + X] using Eq. (|72D: 

^\A^B + D*^D\p)] = T[A^B + X] , ^^^^ „ . (74) 
dM^ ^^'^ ^ ^ 27r|7| |7 + ip|2 ^ ' 

The coefficient of r[A B + X] agrees with Ref. i^]. The corresponding coefficient in 
Refs. fis'] and [i^ is larger by a factor of 2. The origin of this discrepancy is an error by a 
factor of -\/2 in the coalescence amplitude A[D*^D'^ X\ in Refs. |l8| and 21 1. They used 



the universal prediction for this amplitude that was derived in Ref. [l6[. The coalescence 
amplitude is determined by the residue of the pole in the energy for the amplitude for 

^m*o^o ^ ^*o^oi ^ ^ -l/(2Mi2a2). (75) 

^ J 2Mx[^i2 + l/(2Mi2a2)] / 1 12 ; i j 

The universal prediction for this amplitude was first derived in Ref. Il6l and that result was 
used in Refs. [l8| and The error in A\D*^I)'^ — >■ X] in Ref. [l^ came from an error 
in the amplitude J\\D*^I)^ —>■ D*^D^], which was larger by a factor of 2 than the correct 
expression in Eq. ()62|) . The same error in the amplitude A[D*^D^ D*^D^] appears in 
Ref. 0. 

We can exploit the fact that the minimal charm meson model is an effective field theory 
for a more fundamental Lorentz-invariant field theory, namely the Standard Model. This 
implies that C^'V**(m) and C^'V**(m) must have Lorentz-invariant expressions in terms of 
the 4- momenta and polarization 4- vectors of the particles in A and B, the 4- vector Pij_2, 
and the polarization 4- vector whose 3- vector part reduces to e*(m) in the frame where 

Pi% = (Ml+2,0). 

We proceed to deduce the constraints of Lorentz invariance on the discovery production 
process B^ + X and the corresponding production processes for D*^D^ and D^D*^. 

The Lorentz invariance of the fundamental theory is a particularly powerful constraint in this 
case. The only 4-vectors that the short-distance factors can depend on are the 4-momenta 
P^, P^, and Pi^2^ which satisfy Pj^ = P^ + Pi^2^ ^^e polarization 4-vector e^, which 
satisfies P1+2 -6 = 0. The only independent Lorentz scalar that is linear in e* is Pb ■ s*. 
Inner products of the 4-momenta can be expressed in terms of the masses Mb, Mk, and 
Mi+2- Thus Lorentz invariance implies that the T-matrix elements in (j71|) are determined 
by two complex constants C_|_ and C_ defined by 

^+ -e*\m) = C+PB-s*{m), (76a) 



2Ao 

^K+,i _ pK+, 



(m) = C_Pb ■e*(m). (76b) 



2 

The T-matrix elements in Eqs. (j7ip reduce to 

T[B+ ^ K+ + D*^D%p,m)] = ^m,M2 [~C+ A{E,M) + C-] Pb ■ e*{m), (77a) 



T[B+ ^K^ + D'^D^^p^m)] = V4MiM2 ^(Ecm(p)) - C_] • £*(m), (77b) 

T[B+ ^K+ + X(m)] = v/2Mx C+^^P^ ■ e*{m). (77c) 

M12 
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The decay rate for + X can be expressed in the factored form in Eq. ()72j) with 

the short-distance factor 



B+ 



647rM|MiM2 



I a 



(78) 



The invariant mass distributions for the charm mesons in the decays of into +D*'^D'^ 
and + D^D*^ can be expressed in factored forms analogous to Eq. 



dT 
dM 

dV 
dM 



[B+ ^K^ + D'^D^^p)] = rf+p 



7re + ^(7im+P) 
1 



7re + «(7im +P) 



(79a) 
(79b) 



where 



(Mi2/27r)C_/C+ is a complex constant. The nonresonant amplitude 
is required by the renormalization of the operator product expansion. It is completely 
determined by the short- distance coefficients C_ and C+ which are insensitive to the small 
momentum scale |7|. The smallest momentum scale to which they are sensitive is the pion 
mass m^. Since has dimensions of length, the natural order of magnitude of Cg+ is 
l/m^r. Thus if p is of order I7I, the nonresonant terms ±c^+ in Eqs. (fTUI) are suppressed 
by a factor of \'-y\/mT^ compared to the resonant terms. If we keep only the leading terms 
in the expansion in |7|/m, we can set c;^+ = 0. Thus Eqs. (f?^ reduce to Eqs. (f7S|) . In a 
systematic expansion in powers of \^\/m.^, the nonresonant amplitudes ±c^+ in Eqs. fl7^ 
would be retained only if all other effects of the same order in |7|/m7r were also included. 
The effective field theory should be improved to take into account the effective range r_|_ in 
the C = + channel and the scattering length a_ in the C = — channel. One should also 
include terms in the operator product expansion with the operators V^D^D — D^V^D and 
V^DD^ — DV^DK These terms will have operator matrix elements with factors of p^e^ and 
Wilson coefficients suppressed by 

We now apply the factorization formulas to the production process B^ K*^ + X 
and the corresponding production processes for D*^D^ and D^D*^. The only 4-vectors that 
the short-distance factors can depend on are the 4-momenta P^, P^*, and Pi^2 
polarization 4- vector e^, of the K*. Inner products of the 4-momenta can be expressed in 
terms of the masses Mb, M^*, and Mi+2- The independent Lorentz scalars that are linear 
in in e*j^e'' are (Pg ■ £|^)(^b ■ £*), {e*^^, ■ e*), and e^^^/jP^^Pf+a^x^''- Thus the constraint 
of Lorentz invariance reduces the T-matrix elements in Eq. (j7ip to six complex constants 
D±, E±, and F± defined by 



C 



B+ 



2An 



C 



B+ 



c 



e m 



D+Mlg^, + E+Pb^Pb, + P+e^.^^gP^ A+2 ) ^r/<^£*"(m),(80a) 



D-M'^g^u + E^Pb^Pbu + F^efj_yapPBPi+2 



(80b) 



The decay rate for B^ 



K*^ _l_ pg^j^ expressed in the factored form in Eq. (|72p 
The invariant mass distributions for the charm mesons in 
+ D^D*^ can be expressed in factored forms 
analogous to Eqs. ()79|) but considerably more complicated. If we keep only the leading terms 
in the expansions in \'j\/m, these expressions reduce to Eqs. (f73|l. 



with short- distance factor 
the decays of P+ into K*+ + P)*°P)° and K 
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C. Short-distance decay of X 



We now consider the decay of X into a short-distance decay mode C . Examples of such 
decay modes are the discovery mode J/iI)'k^t[~ and J/'ip'^- The expression for the T-matrix 
element analogous to Eq. is 

T[X -^C] = v/2Mx (C^'* {^D'D{Q)\X) + C^'^ {^DD\{])\X)) . (81) 

The operator matrix elements are given by expressions analogous to those on the right 
sides of Eqs. (fTOj) except that e**(m) is replaced by e*(m). The factored expression for the 
T-matrix element is 



r[X(m) - C] = V2Mx ' ll^^ e\m). (82) 

The T-matrix element does not depend on the ultraviolet cutoff if (C"'* + C'"'*)/Ao does 
not depend on A. If we were to consider the process D*^D^ C, we would find that 
_ qc,i j]2ust also be independent of A. Thus the ultraviolet divergences can be removed 
by renormalizations of the Wilson coefficients C*"'* and C*"'*. The decay rate for X C can 
be expressed in a factored form analogous to Eq. 



nX^C] = T^^\^i + ^l\'/\ (83) 

where T'^ is a short-distance factor with dimension of mass. We have chosen the long- 
distance factor to be the same as in Eq. (|48p . The factorization formula in Eq. (|83|) was first 
derived in Ref. 2^ 1. 



In Ref. [2^, the decay rates of X(3872) into J/V'tt+tt", J/^j Tr^ir-ir^ , and J/tpir^'y 

were calculated under the assumption that the decays are dominated by a direct coupling of 
X to J/tp and the vector mesons p and uj followed by the decay of the virtual vector mesons 
into pions and photons. The decay rates were calculated in terms of coupling constants 
Gxtpp and Gx-ipoj and other parameters that were determined by vector meson decays. The 
long-distance scale I7I enters the decay rates only through a factor of 7^/^ in the coupling 
constants Gxi^p and Gx^uj- Thus the decay rates in Ref. |j2J| satisfy the factorization formula 
in Eq. (jHSl). 

The Lorentz invariance of the more fundamental theory provides constraints on the T- 
matrix elements for specific short- distance decay processes. For the decay X — > J/V'7, the 
only 4-vectors the short- distance factors can depend on are the 4-momenta Pi^2 and or 
Pt^ and the polarization 4-vectors and ef^ of the J/ip and the photon. The coefficient 
of elpe*^e'^ must be a 3-index Lorentz tensor. There are 6 independent tensors that can 
be constructed from the 4-momenta, the metric, and the Levi-Civita tensor. Thus Lorentz 
invariance constrains the T-matrix element to be a linear combination of these 6 terms with 



constant coefficients. In the model of Ref. [2^, the assumptions of the direct coupling of 



X to J /ip p and J/ipuj and the vector-meson dominance of the coupling of the photon to 
hadrons were used to reduce the T-matrix to a single term proportional to eo,pupP!^elpe*^ e'^ . 
For the process X —>■ J/ip Tr~^TT~ , the only 4-vectors the short-distance factors can depend on 
are the 4-momenta Pi^2^ -^^5 and Qf^ = P^+ — and the polarization 4- vector e'^ of the 
J I'll). The coefficient of s*^e'^ must be a 2-index Lorentz tensor. There are 8 independent 
tensors that can be constructed from the 4-momenta, the metric, and the Levi-Civita tensor. 
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Thus Lorentz invariance constrains the T-matrix element to be a hnear combination of 
these 8 terms with coefficients that are functions of the two independent Lorentz scalars 
Qtt ■ Pip and Pi+2 ■ P-ip- In the model of Ref. 0|, the assumption of a direct coupling of 
the X to J /tp p was used to reduce the T-matrix element to a single term proportional to 

D. Line shape of X in a short-distance decay mode 

We now consider the line shape of X in the process A ^ B + C , where C is a short- 
distance decay mode of X. An example of such a process is the discovery process for the X: 
+ J/il)Ti^Tt~ . The expression for the T-matrix element analogous to Eq. (|3n|) is 



T[A^ B + C] = Cf^ + {C^'X^'' + C^''C^'')iLo{M, 



c) 

2 



.(C-.C-)(C-.C-)t^^|^. (84) 

The factorized expression for the T-matrix element analogous to Eq. (j3T|) is 

T[A^B + C] = -^ -^^A{Mc) + [C^ + ^C^) 

+ '^{Cl'' - C^'') (C^-^ - C^'') [LoiMc) - Lo(Mi+2)], (85) 

where AC^'^ is 

After renormalization of the coupling constant, A{Mc) is given by the expression in 
which does not depend on the ultraviolet cutoff A. The term Lq{Mc) — Lo{MiJ^2), which 
is given in Eq. (j^ . also does not depend on A. The combinations (C^'* + C^'*)Ao, (C*-^'* + 
C'"'*)/Ao; C^'* ~ and C*"'* — C*"'* cannot depend on A, because they appear as short- 
distance factors in other T-matrix elements such as those in Eqs. (I71|) and (j82j) . Thus the 
T-matrix element in Eq. ()85j) will not depend on A if C^' -|- AC^' does not depend on A. 
Equivalently, the dependence on A can be removed by an additive renormalization of the 
Wilson coefficient C^'*". 

The expression for the rate that follows from the T-matrix element in Eq. (j85|) is very 
complicated and depends on the types of particles in B. The expression simplifies if we 
keep only the leading term in the expansion in |7|/m,r- If \pc\ is of order I7I, the resonant 
amplitude A{Mc) has a factor of order I/I7I while the term Lq{Mc) — Lo(Mi+2) has a factor 
of order |7|. Additional factors of I7I must be accompanied by additional factors of m^r in 
the short- distance factors. Thus the second and third terms on the right side of Eq. ()85|) 
are suppressed by one and two powers of |7|/m^, respectively. If we keep only the leading 
term in |7|/m^, the invariant mass distribution for C in the decay of a single particle A into 
B + C can be expressed in the simple factored form 



(7rc - 

27r 



27r 

7re+ (7im+Pc)^ 



= r^'^ , , ^ Mc > Mi+2. (87b) 
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The invariant mass distribution is continuous at Mc = M1+2 = iVf^io + M^i.o. The factor- 
ization formula in Eq. (jHTj) was first derived in Ref. 23|. In contrast to the corresponding 
factorization formula in the scalar meson model which is given in Eq. ()59|). the short- distance 
factor r^'*^ is not simply the product of the short- distance factors and in Eqs. (jH^ . 
fj73|) . and (jHHj) . The reason for this is that the short- distance factors associated with the 
initial and final states in the T-matrix element in (j85p are connected by the vector index i. 



VIII. SUMMARY 

The X(3872) seems to be a hadronic molecule consisting of a C = + superposition of 
D*^D^ and D^D*^ that are weakly bound in the S-wave channel. The binding energy and the 
width of the X can be conveniently expressed in terms of the complex binding momentum 
7 defined in Eq. (jlUj) . The smallness of 7 compared to the natural scale m^r together with 
the S-wave nature of the bound state imply that the X has universal properties that are 
completely determined by 7. The separation of scales between I7I and m,r can be exploited 
through factorization formulas for the production and short-distance decay rates of X. The 
factorization formulas express these rates as the sum of products of short- distance factors 
that are insensitive to 7 and long-distance factors that are completely determined by 7. 

We have shown how the factorization formulas can be derived using the operator product 
expansion for a low-energy effective field theory for the charm mesons, such as the minimal 
charm meson model. Using the operator product expansion, the rates are expressed as sums 
of products of Wilson coefficients and matrix elements of operators in the effective field 
theory. In the minimal charm meson model, the matrix elements can be calculated nonper- 
turbatively and they depend on the ultraviolet cutoff A. Some of the dependence on A can 
be removed by the renormalization of the coupling constant. This can be accomplished con- 
veniently by eliminating the bare coupling constant in favor of 7. The remaining dependence 
on A can be removed by renormalization of the Wilson coefficients in the operator product 
expansion. After eliminating all dependence on A, the rate can be expanded in powers of 
|7|/m,r- The leading terms in the expansions are very simple. The leading terms in the rates 
for the production processes A^B-\-X,A^B + D*^D^, and A ^ B -\- D^D*° are given 
in Eqs. (j72p and (j73|) . The leading term in the rate for the short distance decay process 
X ^ C is given in Eq. (jHHj) . The leading term for the line shape of X in the short- distance 
decay mode C is given in Eq. (jHZj). 

Our derivation of the factorization formulas using the operator product expansion makes 
it clear how these leading order results can be extended systematically to higher orders in 
|7|/m7r. This requires improving the effective field theory and including higher dimension 
operators in the operator product expansion. If accuracy to n^^ order in |7|/m^ is desired, 
the effective field theory must describe the scattering of charm mesons to that accuracy 
and the operator product expansion must include operators with up to n gradients. As the 
order in \'y\/m.^ increases, there are an increasing number of parameters in the effective field 
theory and an increasing number of Wilson coefficients in the relevant terms of the operator 
product expansion. The difficulty of determining all these parameters phenomenologically 
may limit the utility of the expansion to low orders in j/mTr- 

Our derivation of the simple factorization formulas in Eq s. (1721) . (1731). () 83 p . and ()87p is 
conceptually cleaner than the previous derivations in Refs. [l8|, [2l|, and [2^. Those previous 
derivations were awkward in that they required taking the limit A — >■ cxo while also exploiting 
the fact that the natural scale of the ultraviolet cutoff is m.^. In the present derivation, all 
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dependence on A is removed analytically through renormalization of the coupling constants 
and through renormalization of the Wilson coefficients in the operator product expansion 
without taking the limit A — oo. In a subsequent conceptually independent step, the rates 
are expanded in powers of |7|/m7r. The leading terms in this expansion give the simple 
factorization formulas. 
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